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O . ABSTRACT 

O 

We consider in the unidirectional approximation the propagation of an ul- 
tra short electromagnetic pulse in a resonant medium consisting of molecules 

characterized by a transition operator with both diagonal and non-diagonal 

C ■ 

^ ! matrix elements. We find the zero-curvature representation of the reduced 

r> I Maxwell-Bloch equations in the sharp line limit. This can be used to develop 

the inverse scattering transform method to solve these equations. Finally 
we obtain two types of exact traveling pulse solutions, one with the usual 
exponential decay and another with an algebraic decay. 
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1 Introduction 

If the duration of an electromagnetic pulse is less than all the relaxation times 
in a medium, then the propagation of such a pulse is accompanied only by 
stimulated absorption and re-emission. These electromagnetic pulses we will 
refer to as ultra-short (USP). Usually the study of the propagation of USP 
within the framework of two-level atoms [|l], 0] assumes that the diagonal 
matrix elements of the dipole moment operator are zero (see also the reviews 
[^, Q ^). This model could also describe the very short electromagnetic 
pulses containing a few optical cycles up to a half cycle P, |^, ||, [^ . 

An important approximation for the electromagnetic field called the slowly 
varying envelope approximation, assumes that the pulse envelope is a function 
slowly varying in space and time in comparison to the carrier monochromatic 
wave. In this approximation the Bloch equations describing the evolution of 
two-level atoms and the equation for the pulse envelope can be solved by the 
inverse scattering transform (1ST) method [0, ^ |12]. There are soliton 
and multi-solitons solutions, which describe the propagation and (elastic) 
interaction between USP in a resonant medium. 

Second in importance is the unidirectional wave approximation. This cor- 
responds to the propagation of the electromagnetic wave only in one of the 
two possible directions. The Maxwell and Bloch equations in this approxi- 



mation allow to consider the USP propagation without imposing limitations 
on the pulse duration. As in the first case, the reduced Maxwell-Bloch equa- 
tions can be solved via the 1ST method. Here the soliton corresponds to 
an extremely short pulse of the electromagnetic field, which contains half 
an optical cycle. The breather solutions correspond to pulses containing a 
few optical cycles. As already noted, two-level atoms media have been con- 
sidered with zero diagonal matrix elements of the dipole moment operator. 
In a medium with a linear Stark shift, this assumption should be revised. 
The constant diagonal matrix elements of the dipole moment operator can 
be also induced by the external constant electric field. In an electromagnetic 
wave the resonant levels of this medium will be shifted. In analogy to a Kerr 
medium, for which a high frequency electromagnetic field results into a high 
frequency Kerr effect, a medium possessing a high frequency Stark effect may 



be referred to as a Stark medium. In Ref. [|13[ the interaction of a few-cycle 
electromagnetic pulses with a two-level Stark medium was considered in de- 
tail. Unlike |T3| here we consider a simpler model, where we assume that 



all the matrix elements of the dipole moment operator are parallel to the 
linearly polarized electric field vector, neglect relaxation processes and the 
inhomogeneous broadening of the resonant lines, and assume unidirectional 
waves. This specific model differs from the standard two-level atoms model 
in a minimal degree. This results in a new kind of steady state propagation 
of the USP of extremely short duration. Furthermore, this model admits a 
zero-curvature representation, which can be used as a base to obtain exact 
solution by the 1ST method. After describing the reduced Maxwell-Bloch 
equations in section 2, we give a zero-curvature representation in section 3 



and find steady state pulse solutions in section 4. 



2 The model and the principal equations 



Following ||rT| we consider a plane electromagnetic wave propagating into a 
resonant medium consisting of molecules characterized by the operator of the 
dipole transition between resonant energy levels and let this operator have 
both non-diagonal and diagonal matrix elements. In the two-level approxi- 
mation the Hamiltonian of the considered model can be written as 

/ 



H = — 



-1 \ / diiE duE 

y 1 y (^ d2lE d22E 

where E is the strength of electric field of the electromagnetic wave. The 
polarization of the medium is P = uaP-, where ua is the density of the 
molecules and the polarizability p is given by the expression p =tipd = 
Pudu + P22'^22 + Pi2'^2i + P2i'^i2- Wc will cousidcr the propagation of short 
electromagnetic pulses, for which all relaxation processes can be neglected. 
Hence, the density matrix p obeys the constraint pn + P22 = 1- Taking this 
relation into account, the polarizability can be written as 

1 1 

The evolution of the elements of p is given by the equation ihdp/dt = Hp — 
pH and yields the Bloch equations. 

In a scalar form the Maxwell equations lead to 

d'^E 1 d'^E Attua d^ 



dz^ c^ dt^ (? df^ 



(P) , (1) 



where the angular brackets stand for averaging over all molecules. 

It is convenient to introduce the Bloch vector by defining its components 

ri = Pi2 + p2i^ ^2 = -iipu- P2i)^ r3 = p22-Pu- 

On can choose a constant phase for the matrix elements of the density matrix 
p and the dipole operator d so that du = '^21 • The polarizability p of the 
resonant medium can be written in terms of the Bloch vector components as 
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P = - (du + d22) + - id22 - dii) ra + di2ri. (2) 

In this expression the first term corresponds to the constant polarizability of 
the molecules. As it gives no contribution to the radiation, it may be omitted. 
Furthermore, after averaging over all atoms it must be {du + ^22) = if one 
assumes that there is no constant polarization of the medium in the absence 
of an electromagnetic field. 

So the total system of equations describing the model under consideration 
can be written in the following form 

d^E 1 d^E AnuA d^ /I , ^ 

- [d22 - du) rg + di2ri 



Or I 

— = -[ujo+ {dii - ^22) E/h] r2, 

^ = - [^0 + (^11 - d22) E/h] n + 2{dE/h)r3, 

^ = -2{dE/h)r2. 

This system differs from the well-known Maxwell-Bloch equations 0- by 
terms containing the parameter {du — ^22)- If ^a \dii — ^22! and ua 1^12! are 
small, then one can neglect the refiected wave and consider the unidirectional 



propagation of the electromagnetic wave |^, ^]. The system of the Maxwell- 
Bloch equations is then reduced to 

dE IdE 27rnAd /I , ^ ^ , \ .ox 

-^ = -[ujo + idu-d22)E/h]r,, (4) 

^ = - [cuo + (rfii - d22) E/h] n + 2idE/h)rs, (5) 

^ = -2(rfE/;^)r2. (6) 

Let us consider here only the sharp line limit of the model. Introduce new 
dimensionless variables and field 

T = uJo{t- z/c), C = z/Lab, q = 2dE/huJo, 

where L~i^ = 4:7mAdl2{fic)~^ and the parameter /i = {dn — ^22) /2(ii2- The 
reduced Maxwell- Bloch equations (2.3)-(2.6) take the form: 

ri,r = -(1 + /ig)r2, r2,r = (1 + /ig)n + qrs, rs^r = -qr2, (7) 

q,C = -(ri - /irs),^, (8) 

The last equation can be rewritten as 

q,C = r2- (9) 

We will assume the general initial conditions where the medium is initially 
at rest so that q = 0, ri = r2 = 0, r^ = —1, at r — > —00. 

From the Bloch equations (^ we obtain ('"i + '"i + ^3) t ~ 0- ^^^ using 
the initial conditions we get the following integral of motion 

rl + rl + rl = 1. (10) 



We now proceed to give the zero- curvature representation of the system (|^- 
d^). We will see that the /xg term which makes the difference from the classi- 
cal reduced Maxwell-Bloch equations [^ |^, |1^ will generate complications of 
traditional 1ST method. We will also find two types of traveling pulse solu- 
tions which differ from the ones of the two classical reduced Maxwell-Bloch 
systems. 

3 The zero- curvature representation 

The inverse scattering transform method is based on considering the pair of 
linear equations |T5| : 



^,. = Uij, V^,^ = V^, (11) 
where 

The compatibility condition for these equations is t/^ = V^^ + VU — UV, 
which yields the following system of equations 

A,, + U21B - UuC = f/n,c, (12) 

B,r + iU22 - Un) B + 2U^2A = Uu,o (13) 

C,r - {U22 - f/ll) C - 2U21A = U2U- (14) 

It should be pointed that the reduced Maxwell-Bloch equations (|3)-(i|) can 
be represented in another form: 

iiP,^ = -(1/2)(1 + ^^q)i,^ - (l/2)g^2, (15) 



^^2,. = -(l/2)gV^i + (1/2)(1 - /ig)^i, (16) 

9,C = ^ (^1^2 - V'iV'2) ' (17) 

where the components of the Bloch vector are constructed from the functions 
ijj according to the rule: 

^1 = (V'tV'2 + V'iV'2) ) ^2 = ^ (^IV'2 - V'iV'2) > ^3 = \1jJ2f - lipif . 

This suggests that the U-matrix of the zero-curvature representation for the 
reduced Maxwell- Bloch equations (0)-® could be written as 

\^ 9i(A), a - /2(A), ) 

where /i,2(A) and (71, 2(A) are unknown functions of the spectral parameter 
A. Taking into account this U-matrix, the compatibility conditions ([T^)-([l^ 
can be written as 



A,. + g,i\)qB - g,{\)qC = /i(A)g,c, (19) 

E. + (2r/ - f{X)q) B + 2g^{X)qA = (?i(A)g,o (20) 

C, - {2^ - f{X)q) C - 2g^{X)qA = g2{X)q,c, (21) 

where /(A) = /i(A) -|- /2(A), and t] = iX. The system of equations ([T9|)-(pTD 



must coincide with the reduced Maxwell-Bloch equations (0)-(|). 

It should be noted that the system of equations (|3)-(iD contains the partial 
derivatives only of the Bloch vector components with respect to r. Hence 
one can find B and C from (^) and ( pT]) using the following expansions 

B = bin + &2^2 + hrs, C = CiVi + C2r2 + c^r^, 



where 61 2,3 and Ci 2,3 are unknown functions of A. Substituting these expan- 
sions into equations ( pO|) and (pT]), and equating the coefficients of ri,r2, r-^ 
and q in the resuhing equations, we obtain the following relations between 
the 61,2,3 and Ci,2,3 



62 + 2r/6i = 0, -hi + 2r]h2 = gi, 63 = 0, 

C2 - 2r]Ci = 0, -ci - 2r]c = 5-2, C3 = 0, 
2giA + hi {-fir2 - fr^) + 62 (yuri + rg - /r2) = 0, 
- 2g2A + ci (-/ir2 + fn) + C2 (/xri + rg + /r2) = 0. 
From (|22|) and (El) it follows that 



(22) 

(23) 
(24) 
(25) 



hi 



-9i 



2r/^i 



Cl 



-^2 



C2 



-2?75f2 



1 + 4^2 ' 1 + 4r/2 ' 1 + 4r/2 ' 1 + Ar]^' 

Substitution of these expressions into (^) and ( [25| ) leads to two equations 
for the matrix element A. From (3.14) one gets 



-1 



2A + (1 + 4r7"j [2r]rs + (2r//i + /) n + (/i - 27^/) ra] = 0, 

and from (3.15) one gets 

2A + (1 + 47^2)"' [2r]r; + {2r]fi + f)r,-{^i- 2r]f) r^] = 0. 

Both expressions agree if /i — 2r]f = 0. Thus the auxiliary function /(A) is 
specified: /(A) = /x/2?7, and the expressions for B, C and A are determined: 



B 



91 



A 



(-n + 2r/r2) , C 
-1 



^2 



l + Ar]^ 1 + 47^2 



(-ri - 2r7r2) 



1 + 4^2 



r^rg + /i I r] + — I ri 



(26) 
(27) 
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To determine the auxiliary function (71, 2(A), equation (3.9) has to be used. 
We assume the equahty /i(A) = /2(A) so that /i(A) = yu/4?7. Substitution of 
the equations (|26D and p7| ) into (|19D leads to the following relation 

Vr3,T + /i I r] + — j n^r + — (1 + 4?]^) g,^ = gig2^qr2- 

Using equations (^ and @) this becomes 

r/ + ^ (1 + 4r/2) + Ar]g,g2 = 0. (28) 

Now, we can require any correlation between gi and g2i for instance, gi{\) = 



—(72(A) = g{X)- From ( pSj) it then follows immediately that 

nl/2 



.(A) - \ 



1 + 7^(1 + 4^^ 



(29) 



Some times it may be more convenient to take gi{\) = (72 (A) = «(7(A). 

Thus, we have the U-V-matrices for a zero- curvature representation of 
the reduced Maxwell- Bloch equations (0) and (^ describing the USP prop- 
agation in Stark medium assuming a model of two-level atoms. They are 



/ 



U 



-7] + {^i/Arf) q 
I + {^ / 27] f {I + Arf 



1/2 



1 + (^/2r7)^ (1 + At,'] 



1/2 



(30) 



V 



2(1 + Ari^ 



■X 



-2 [r]r3 + fi{r] + (1/4?])) n] 
1 + {fx/2rif (1 + Ari^)] ^'^ {2r]r2 + n 



1 + in/2rif (1 + Ari^)] ^''^ {2i]r2 - ri) 



2 [r/rs -h /i (r/ + (/i/4?7))ri] 



(31) 
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where t] = iX. It should be noted that if fi — > 0, then we obtain the well- 
known U-V representation for the classical reduced Maxwell- Bloch equations 

(see 



O'Hl)- 



4 Steady state solutions of the Maxwell-Bloch 
equations 

To obtain the equations describing the propagation of stationary USP one 
should assume that the components of the Bloch vector and the normalized 
pulse envelope depend on only one variable t] = t — (]( = ujQ{t — z/V). 
Under this assumption the system (|^ and (|^) transforms into the system of 
ordinary differential equations 

n^rj = -(1 + /ig)r2, r2,^ = (1 + Aig)n + ^'"a, '"s.r, = -qr2, (32) 

(3q,v = -r2- (33) 

Let us consider the solution of equations (|3^) and (^3|) describing a soli- 
tary steady state wave with the boundary conditions q = 0,ri = r2 = 0,r^ = 
—1 at \ri\ — > oo . 
From (|32|) it follows that 

r2 = -(ri -yurs),^ (34) 

We now take (P3| ) and the boundary conditions into account to get 

n = Pq + fx{l + rs). (35) 

From the last equation of (^) and equation (p3D we obtain 

r3 = -l + /?gV2. (36) 

11 



We define q^ = 4//5 to renormalize q as w = q/qo and write all components 
of the Bloch vector as ri = iw/qo + 2fiw, r2 = —iw^rj/qo, '"s H — 1 + 2w^. The 
substitution of the Bloch vector components into expression (|10[) leads at 
once to the equation for w: 



' dw 



2 



-v2 \ ^2 



l\w'- fiqow' - f (l + f^') w'- (37) 



\drj ) \4 / 4 

This equation has a solution in the form of a solitary wave if the factor in 
the first term is positive. Let us denote it as go/4 — 1 = 0~'^ and introduce 
the new variable 



,2 



y = 2 [qQW\Jl + li 
In terms of this variable equation (^) takes the form 

Let be M = y — fiO'^ /{I + /i^)^/^.. Then equation (|38|) can be rewritten as 

du\ / 2 d2\i/2 



"U'-'--^ 



where 

52 = ^M 1 



2 /1 , At ^ 



The solution of this equation is u{ri) = 5 cosh [{rj — tjq) /6], with rj^ as the 
integrating constant. Taking this result into account one can obtain that 

eUi + /i2(i + e^) cosh [{ri - rio) /e] + fxe\ 

This expression represents a one-parameter family of solutions of equations 
([3^ ) and (|33D , with 6 as parameter. It is convenient to introduce the duration 

12 



of a stationary USP ( PD| ) tg by the relation 6 = tgUJo, so that q (77) may be 
written as 

Q iv) = ( I r- (40) 

1^1 + /^n + {^tsUJ^Y) cosh [it - z/V - to) As] + litsOJo] 

Thus one can obtain the expression for the electric field of this extremely 
short steady state pulse: 

E, (t, z)= ^. (41) 

1^1 + /i^l + (/it.CJo)^) cosh \{t - z/V - to) /ts] + /it.CJoj 

From the definition of r] one can obtain the velocity of propagation V of the 
USP in the laboratory reference frame 

This expression for the velocity of the propagation of stationary pulse is 
identical to the one found in ||^. Hence the linear high frequency Stark effect 
does not influence the rate of propagation of the steady state USP. 



The expression for the normalized electric field strength (|39|) has been 
obtained under the condition that 9^^ = q^/A — 1 > 0. Let us consider the 
case 9 = 0, i.e., go = 2. Equation (|37|) takes the form 

If fi is positive we have only trivial solutions of this equation, w = 0. If /x is 
negative there are nontrivial solutions. We rewrite (^) as 

dr] I 



'^^^^ ' [2 |/i| - (1 + /i>] . (44) 



w^ 
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and introduce -u = 2 |/i| / (1 + /i^) w so that the equation becomes 

The solution of this simple equation is 

Hence 

^ = Ti;^V + ITT? <"-"»'/• 

and 

In terms of physical values, this algebraic steady state USP can be rep- 
resented as 

Eai (t, z) = — 2 , (46) 

where the amplitude of USP Em, its duration tau and velocity Vai are defined 

as 

_ 4:huJo{d22 - du) _ C^o(rf22 - c^ii) JL - ^4 ^'^^Ad'^ 
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^MJ2 + {d22-diif " ^MJ2 + {d22-dnf ^-i ^^o 

(47) 

The possibility of the propagation of such an algebraic solitary wave is the 

distinguishing feature of the Stark medium under consideration here. 



5 Conclusion 

We have introduced and analyzed a model for the propagation of ultra-short 
electromagnetic pulses moving in one direction in a Stark medium. This is 

14 



described by two-level atoms taking into account the high frequency linear 
Stark shift of the energy levels. Two families of exact analytical solutions of 
the reduced Maxwell-Bloch equations have been found, which correspond to 
steady state pulses propagating with different pulse widths. We find algebraic 
solitary waves apart from the usual exponentially decreasing ones. 

It was found that the system of reduced Maxwell-Bloch equations admits 
a zero-curvature representation in the sharp line limit. Unlike the well-known 
Zakharov-Shabat spectral problem the one obtained here shows complicated 
analytical properties. Nevertheless we assume that it can be used to develop 
the inverse scattering transform method to solve these equations. 

There is another class of stationary waves - cnoidal waves amongst the 
solutions of the Maxwell-Bloch equations (0)-(||), or (0) and (|]). These are 
periodic continuous waves, different from the solitary waves described above. 
Since equations (0) and (P) are valid when the duration of the wave is less or 
much less than the relaxation times of the atomic subsystem, cnoidal waves 
seem to be mathematical objects lying beyond the physical meaning of the 
original equations. 

6 Acknowledgment 

One of the authors (A.I.M.) is grateful to the Laboratoire de Mathematiques, 
INS A de Rouen for hospitality and support. 



References 



15 



[1] S.L. McCall, E.L. Hahn, Phys. Rev. Letts. 18 (1967) 908 ; Phys. Rev. 
183 (1969) 457. 

[2] L. Allen, J.H. Eberly, Optical Resonance and Two-Level Atoms (Wiley, 

N.Y. 1975). 

[3] G.L.Lamb, Jr., Rev.Mod Phys. 43, (1971) 99. 

[4] R.K. Bullough, P.M. Jack, P.W. Kitchenside, R. Saudders , Phys.Scr. 
20 (1979) 364. 

[5] A.I. Maimistov, A.M. Basharov, S.O. Elyutin, Yu.M. Sklyarov, Phys. 
Rept. 191 (1990) 1. 

[6] R.K. Bullough, F. Ahmad, Phys.Rev.Lett. 27 (1971) 330. 

[7] A. E. Kaplan, S.F. Straub and P. L. Shkolnikov, J.Opt.Soc.Amer. B 14 
(1997) 3013. 

[8] A.I. Maimistov, Quantum Electron. 27 (1997) 935. 

[9] A.I. Maimistov, Quantum Electron. 30 (2000) 287. 

[10] G.L.Lamb,Jr. Phys.Rev. A, 9 (1974) 422. 

[11] M.J.Ablowitz, D.J.Kaup, A.C.Newell, H.Segur, J.Math.Phys. 15 (1974) 
1852. 

[12] D.J. Kaup, Phys.Rev.A, 16 (1975). 

[13] L.W. Casperson, Phys. Rev. A 57 (1998) 609. 



16 



[14] J.L. Eilbeck, J.D. Gibbon, P.J. Caudrey, R.K. Bullough, J. Phys. A 6 

(1973) 1337. 

[15] M.J.Ablowitz, D.J.Kaup, A.C.Newell, H.Segur, Stud.Appl.Math. 53, 

(1974) 249. 



17 



